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Abstract 

The production of W + W~ pairs is analysed in a data sample collected by ALEPH at a 
mean centre-of-mass energy of 182.7 GeV, corresponding to an integrated luminosity of 
57 pb _1 . Cross sections are given for different topologies of W decays into leptons or 
hadrons. Under Standard Model assumptions for the W-pair production and decay, the 
W-pair cross section is measured to be 15.57 ± 0.62(stat.) ± 0.29(syst.) pb. Using also 
the W-pair data samples collected by ALEPH at lower centre-of-mass energies, the decay 
branching ratio of the W boson into hadrons is measured to be B(W — > hadrons) = 68.93± 
1.21(stat.) ± 0.51(syst.)%, allowing a determination of the CKM matrix element |V^ S | = 
1.043 ±0.058(stat.) ±0.026(syst.). The agreement of the cross sections with the Standard 
Model prediction allows a limit to be set on the W decay rate to undetectable final states. 
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1 Introduction 



This letter presents results on W-pair production using data collected with the ALEPH 
detector at centre-of-mass energies ranging from 180.8 to 183.8 GeV, during the 1997 data 
taking period. 

The experimental conditions and data analysis follow closely those used in the cross 
section measurement at lower LEP2 energies. As they are already described in detail 
in [|lj, attention will be focused on changes in selection procedures other than a simple 
rescaling of cuts with the increased collision energy. 

At these energies the cross section has a much smaller dependence on the W mass with 
respect to the threshold productions, and the total production rate constitutes a test of 
the Standard Model (SM). 

About 900 W pairs are expected in the collected data sample, allowing improved direct 
determinations of the W hadronic and leptonic branching ratios. The hadronic branching 
ratio is sensitive to the yet poorly known coupling |V^ S | of the W to cs pairs. 

A detailed description of the ALEPH detector can be found in Ref. and of 
its performance in Ref. ||. The luminosity is measured with small-angle Bhabha 
events, using lead-proportional wire sampling calorimeters [|J, with an accepted 
Bhabha cross section of approximately 4.6 nb |J. An integrated luminosity of 
56.81 ± 0.11 (stat.) ± 0.29 (syst.) pb _1 was recorded, at a mean centre-of-mass energy of 
182.66 ± 0.05 GeV @. 

2 Physics processes and definition of the W-pair 
cross section 

To lowest order within the Standard Model, three diagrams contribute to W-pair 
production in e + e~ annihilations: the s-channel 7 and Z boson exchange and the t- 
channel z/ e exchange, referred to as CC03 diagrams. Each W is expected to decay 
rapidly into a quark- ant iquark pair (W~ — > ud,cd, us, cs, ub, cb) or a lepton-neutrino pair 
(W~ — > e~V e , fi~V^,T~V T ). Therefore, each CC03 diagram leads to an experimentally 
accessible four-fermion final state. However, many other Standard Model processes can 
lead to the same four-fermion final states as W-pair decays, interfering with the CC03 
diagrams. 

In the following, all signal cross sections are defined by the production of four-fermion 
final states only through two resonating W bosons, and will be referred as CC03 cross 
sections. The effect of the non-CC03 diagrams is corrected for with additive terms, 
obtained by comparing Monte Carlo (MC) simulations including only CC03 processes 
with all Standard Model four-fermion final states compatible with W-pair decays (WW- 
like four-fermion final states), following the same procedure as in Ref. [lj. Such Monte 
Carlo additive terms are referred to in the following as "4f-CC03 corrections". They 
are computed separately for all different WW-like four-fermion final states /1/2/3/4, as 
(e4f<74f — eccosCccoa)^ 2 "^ 4 , where e and a refer to MC selection efficiencies and generator 
cross sections, on samples of /1/2/3/4 final states generated with the full four-fermion 
productions (4f subscript) or with W-pair decays only (CC03 subscript). All the signal 
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efficiencies quoted in the following are to be interpreted as CC03 efficiencies, evaluated 
on such Monte Carlo samples. 

Two Monte Carlo event generators are used to simulate the signal events. Samples 
of events are generated with different W masses, both for CC03 diagrams and for all 
WW-like four-fermion diagrams, with KORALW 0. A comparison sample is generated 
with EXCALIBUR || with and without colour reconnection effects (following the ansatz 
of Ref. H). The KORALW samples serve to determine the four-fermion and CC03 
efficiencies used to obtain the final result. Other KORALW samples are used to check the 
m\v dependence of the selection procedures and of the four-fermion to CC03 correction. 
The EXCALIBUR samples are used as a check of the Monte Carlo simulation of the physics 
processes, and to assess the effects of colour reconnection. 

The PYTHIA 5.7 |jTD| Monte Carlo program is used to generate background events 
coming from qq, ZZ, Zee and Weu processes. Other background events are generated 
with PH0T02 [IT for gamma-gamma interactions, KORALZ [[EJ for fi and r pair production, 
BHWIDE [13] and UNIBAB [14| for Bhabha scattering events. For the fully leptonic selections, 
non- WW-like four-fermion background processes with two visible leptons (as £ 't'vgvg) are 
generated with EXCALIBUR. Care is taken to avoid double counting of backgrounds with 
different generators and to avoid counting as backgrounds those WW-like four-fermion 
processes that are part of the 4f-CC03 correction. Monte Carlo samples corresponding 
to integrated luminosities at least twenty times as large as that of the data are fully 
simulated for all background reactions. 



3 Selection of W-pair candidates 
3.1 WW — ► tvlv events 

The selection of fully leptonic W-pair decays is an update of the two selections used for 
the cross section and branching ratio measurements at 161 and 172 GeV 0, to the 183 
GeV energy regime. 

In the first selection, which does not make use of lepton identification criteria, two 
changes are introduced. The lepton-lepton acoplanarity is required to be smaller than 
175° and the energy of the leading lepton must be smaller than 86 GeV. In the data, 47 
events are selected with this analysis. 

In the second selection lepton identification is used to classify events in six different 
di-lepton channels (ee, e/z, er, /i/x, /ir or rr), according to the flavour of leptons, 
and optimised cuts are applied individually in each channel. The electron and muon 
identification criteria are the same as those used for the leptonic classifications at lower 
energies jl]. A lepton then is classified as a tau either if no lepton identification is fulfilled 
or if the identified lepton has an energy lower than 25 GeV (0.137^). In the data, 57 
events are selected with this second analysis. 

The inclusive combination of the two selections has an efficiency of 71.5% for the fully 
leptonic WW channels, combined assuming lepton universality. For the cross sections 
and the branching ratio measurements, the events selected by the combination of the two 
analyses are classified into six channels according to the di-lepton flavour classification 
used for the second selection. The CC03 efficiencies in the individual iviv channels are 
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given in Table [I] for the inclusive combination of the two selections. The total background 
amounts to 151 fb and is dominated by 77 — > li and non-WW-like four-fermion events, 
as e + e~— > £'£'i/£Vz with £' 7^ £. In the data, the inclusive combination selects 61 events. 

Table 1: Summary of results of the different event selections on Monte Carlo and Data 
events. Efficiencies are given in percent of CC03 processes. The Tvqq selection column 
refers to the performance of the exclusive selection. The qqqq column refers to events 
with a NN output greater than -0.2 . The listed backgrounds do not include the 4f-CC03 
corrections. In the qqqq column the backgrounds include non-qqqq WW decays. 
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Beam related background and detector noise are not fully reproduced by the simulation 
and affect the efficiency of the photon veto cuts. Random trigger events are used to 
compute the inefficiency introduced. The result is (4.3±0.8)% for the low angle photon 
veto and is applied as a correction factor lowering the Monte Carlo efficiencies. The 
systematic error takes into account the uncertainty on the inefficiencies. For the large 
angle photon veto these effects are not corrected for, but a systematic error of ±1.7% is 
assigned to them, again from the study of random trigger events. Lepton identification 
efficiency has a 2% systematic uncertainty, to cover residual discrepancies in identification 
efficiency and purity between LEP1 data and Monte Carlo. The final sources of systematic 
errors are listed Table [2| The effect on the total cross section amounts to 0.04 pb and is 
dominated by the Monte Carlo statistical errors in the subtraction of backgrounds and in 
the signal efficiencies. 

A likelihood fit is applied to determine individual cross sections for each fully leptonic 
decay channel using efficiency matrices for CC03 and backgrounds, as in Table [p. WW-like 
four-fermion processes are taken into account in the 4f-CC03 correction for each individual 
channel, given by (e 4 fO"4f — ecco30"cco3)^ l/ ^ l/ - For all channels together, this correction 
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Table 2: Systematic errors summary in fb. 
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amounts to +6 ± 19 fb, where the uncertainty comes from Monte Carlo statistics. 
The results of the fit, where the cross sections are constrained to be positive, yield 
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± 0.01(syst.) pb, 
±0.01(syst.) pb, 
± 0.03(syst.) pb, 
± 0.02(syst.) pb, 
± 0.03(syst.) pb, 
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where the statistical errors correspond to a variation of the likelihood logarithm of 0.5 
around the fitted maximum. The systematic errors are obtained by varying all input 
parameters in the fit according to their uncertainties. For cr (WW — > tvtv), the systematic 
uncertainty refers to the positive statistical error, rather than to the fitted central value. 

The fully leptonic cross section is obtained with the same fit, assuming lepton 
universality, and yields 



cr(WW -> tvtv) = 1.38 ± 0.20(stat.) ± 0.04(syst.) pb. 



3.2 WW — > £vqq events 

The three selection procedures developed for the lower energy measurements are applied. 
One selection is optimised for WW events with electrons or muons and requires an 
energetic identified electron or muon, while the other two are designed for Tvqq events, 
based on global variables and topological properties of events. 
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3.2.1 WW — > ei/qq and WW — > /xi/qq selections 

The selection of ez/qq and /zz/qq events is basically unchanged with respect to the previous 
analysis fl|]. After the preselection of hadronic events based on charged energy and 
multiplicity, the charged track with the highest momentum component antiparallel to 
the missing momentum is chosen as the lepton candidate. The same electron or muon 
identification criteria as for the fully leptonic channel are required for this lepton candidate 
track, as well as an energy of at least 15 GeV. The probability for an event to come from 
a signal process is determined from the energy and isolation of the lepton and the total 
missing transverse momentum. Selected events are required to have a probability greater 
than 0.55, both for electron and muon candidates (Fig. |I|). 

The CC03 selection efficiencies are 81.5% for the ez/qq and 88.8% for the /iz/qq channel, 
with a total background of 158 fb; 127 electron and 113 muon events are selected in the 
data. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Prob(evqq) Prob((lvqq) 

Figure 1: Probability distribution of preselected events for the ez/qq (a) and the /zz/qq (b) 
selections. The arrows show the cuts above which events are kept in the selection. 

3.2.2 WW — > ri/qq events 

The selection of rz/qq events is based on global event variables and is combined with a 
topological selection designed to identify the r jet. The common preselection is described 
in Ref. [|IJ. In the global analysis the acoplanarity is required to be less than 175°, the 
energy in a wedge of half-angle 30° centred on the missing momentum direction in the 
plane transverse to the beam axis smaller than 12.5%^, the estimated energy of the 
"primary" neutrino smaller than 65 GeV, and the visible mass smaller than 140 GeV/ c 2 . 
In the topological analysis the higher energy of the two quark jets must be smaller than 
65 GeV and the acolinearity of the quark jets more than 110°. 
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In the data 238 events are selected by the combined analyses, of which 152 are already 
in the ez/qq and /zz/qq samples. The tau exclusive selection, i.e. not counting events already 
selected by the electron and muon semileptonic selections, has an efficiency of 54.1% and 
total background of 114 fb, dominated by qq events. 

3.2.3 Results 

The inclusive combined efficiencies for the three semileptonic decay channels are 87.7% for 
the electron channel, 92.3% for the muon channel and 60.5% for the tau channel, giving 
an 80.2% overall efficiency for WW — > iuqq, with a background of 270 fb (Table H). The 
overall 4f-CC03 correction is +35 ± 36 fb. A total of 322 events are selected in the data. 

The systematic uncertainties on the cross sections (Table |2|) are dominated by the 
background normalisation, and by the uncertainty coming from the cut on energy around 
the beam for ruqq events. Since the selections of semileptonic events is basically 
unchanged, all systematic uncertainties on the signal efficiencies and on the backgrounds 
have been evaluated in the same way as for the lower energy measurements |l||. 

A similar fit as for fully leptonic events is used here with the corresponding matrix 
of efficiencies and backgrounds including the overlaps of the different iuqq analyses. The 
partial cross sections are then extracted from such a maximum likelihood fit to the number 
of events in each selection, leading to 

<r(WW -> ez/qq) = 2.52 ± 0.24(stat.) ± 0.06(syst.) pb, 
ct(WW -> fiuqq) = 2.12 ± 0.21(stat.) ± 0.06(syst.) pb, 
cr(WW -> ruqq) = 2.15 ± 0.31(stat.) ± 0.09(syst.) pb, 

where the systematic errors are obtained by varying all input parameters of the fit 
according to their uncertainties. 

The total iuqq cross section is extracted by means of the same fit, under the 
assumption of lepton universality: 

a(WW -> £uqq) = 6.81 ± 0.40(stat.) ± 0.14(syst.) pb. 
3.3 WW — ► qqqq events 

The analysis of WW decays to four jets is updated from Ref. [IJ and consists of a simple 
preselection followed by a fit to the distribution of the output of a neural net (NN). 

For each event, particles are clustered into four jets with the DURHAM algorithm fill . 
A first preselection cut rejects events with a real Z and large undetected initial state 
radiation, requiring that the total longitudinal momentum be smaller than 95% of the 
difference between the visible mass and the mass of the Z. The sphericity must be larger 
than 0.03. The value Y34 of the jet resolution parameter where a four jet becomes a three 
jet event, is required to be greater than 0.001. A further requirement that none of the four 
jets contains more than 95% of electromagnetic energy rejects qq events with a visible 
initial state radiation photon. 
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Figure 2: Comparison of NN output distributions for data and Monte Carlo after the 
preselection. The points are the data and the open histogram the total Monte Carlo 
prediction. The light shaded histogram shows the expected WW — > qqqq contribution. 

At this level of the selection, 1014 events are selected in the data while 995.4 events 
are expected from all Standard Model processes. This preselection has an efficiency of 
97% on CC03 events and a purity of 40%. 

The input variables for the NN are described in Appendix A, and are related to the 
global event properties, the properties of jets, WW kinematics and the b-tag probabilities 
for the four jets. Figure |2] gives the distribution of the NN output values for all 
backgrounds and signal Monte Carlo events compared with the data. 

The number of signal events is extracted by means of a binned maximum likelihood 
fit to the distribution of the NN output for data events. Only the normalisation of the 
MC signal is allowed to vary in the fit. The 4f-CC03 correction is +10 ± 8fb. 

The fit result is 

(j(WW -> qqqq) = 7.35 ± 0.42(stat.) ± 0.25(syst.) pb. 

Systematic errors are summarised in Table 0. The uncertainty on the variables used 
for the event preselection cuts, and for the NN input (preselection and shape variables), 
has been assessed by reweighting the events in the NN output according to bin by 
bin Data/MC differences in the distributions of such variables. Uncertainties with the 
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WW generator have been evaluated comparing KORALW and EXCALIBUR MC samples, 
EXCALIBUR has also been used to evaluate the effects of colour reconnection. Finally, 
for the QCD fragmentation uncertainty, the HERWIG [16] generator was used to produce 
different samples of WW signal and qq background events. 



4 Total cross section 

The total cross section is obtained from a fit assuming the Standard Model branching 
ratios, the only unknown being the total cross section. The fit is applied to all data 
selected as described in the previous sections, and using the matrices of efficiencies and 
backgrounds of the various analyses, yielding 

a ww = 15.57 ± 0.62(stat.) ± 0.29(syst.) pb. 

Figure |3] shows the total cross section measured as a function of the cm. energy. 



5 Branching ratios and V cs 

The same fit as for the total cross section is performed adding together the data samples 
collected at 161, 172 and 183 GeV centre-of-mass energies. 

Without assuming lepton coupling universality, the six unknowns are the three 
individual leptonic branching ratios and the three total cross sections at 161, 172 and 
183 GeV. The hadronic branching ratio is set to 1 — B c — — B T . The fitted leptonic 
branching ratios are 

B(W^ei/) = 11.15 ± 0.85(stat.) ±0.24(syst.)%, 
B(W -»• fiv) = 10.06 ± 0.78(stat.) ±0.21(syst.)%, 
B(W->™) = 9.76 ± 1.01(stat.) ± 0.33(syst.)%, 

and are consistent with lepton universality and the Standard Model expectations. Due 
to cross-contaminations in the identification of W decays to tv versus eu or fiv, the 
measured B(W — > tv) is 26% anticorrelated with B(W — > ev) and 22% anticorrelated 
with B(W -»• fiv). 

If lepton universality is assumed (B e = B^ = B T = (1 — B q )/3), a fit to B(W — > qq) 
and the three total cross sections as unknowns yields the measurement of the hadronic 
branching ratio 

B(W -> qq) = 68.93 ± 1.21(stat.) ± 0.51(syst.)%. 

This result can be expressed in terms of the individual couplings of the W to quark- 
antiquark pairs: 

mX ~ " <IMi (IKal 2 + IKdl 2 + IKsl 2 + IKsl 2 + IKbl 2 + |K b | 2 )(l + a s «)A0 



1 - B(W -> qq) 
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Figure 3: Measurements of the W-pair production cross section at three centre-of-mass 
energies, compared to the Standard Model prediction from GENTLE j^Ej] for the world 
average value of the W mass p^/ . The two curves correspond to the 80 MeV/ c 2 error 
on m\y ; the shaded area around the two curves represents the theoretical error (±2%) 
on the GENTLE calculations. 



Using the world average value of a s (m|) evolved to together with the other measured 
CKM matrix elements 



17 allows a constraint on the least well measured CKM matrix 



element, 



\Vcs\ = 1.043 ± 0.058(stat.) ± 0.026(syst. 
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6 Invisible W width 



The cross section given in the previous sections are based on the assumption that the 
W decays exclusively into the known quarks and leptons. Scenarios of new physics have 
been advocated whereby the W could decay into a charged particle with momentum below 
detectability, P < 200 MeV/ c |20fl . These decays would appear as invisible, in that events 
containing them would not be selected by the previously described selections. A first 
consequence would be the appearance of events with only one visible W decay. A direct 



search for this scenario has been performed by ALEPH and is described elsewhere |2T . 

A second consequence would be a modification of the visible WW cross section as 
follows. The total width T-yv of the W boson would be increased with respect to the 
Standard Model prediction for decays into known fermions, r^Jf, by the corresponding 
invisible partial width, 

p "pvis i "pinvis 

i-w — 1 w ± 1 w • 



The dominant effect would be a decrease of the visible cross section cr^f w according to 



"WW 

where 



(S™)Vww, 



"pvis 

^vis W -y pinvis 

pvis I pinvis 

1 W "T 1 W 

The total cross section cxww has a non- negligible, but small, dependence upon a change 
of the total width, if one assumes that the couplings of the W to known fermions are 
left unmodified. The above formulae have to be corrected for mixed decays, in which 
one W boson decays to a visible channel while the other W boson decay is undetected. 
These events would not be completely missed by the standard analyses, they have an 
efficiency to be selected by the rvqq selection of 22.4% and 8.9% at 172 and 183 GeV, 
respectively. Such events would not be detected by the other selections. This correction 
term is proportional to r^ vls /r^ s and the above selection efficiency; it represents a change 
in the visible cross section at the per mil level. 

The invisible width T^ V1S is extracted by means of a fit to the measured cross sections 
at y / i=161, 172 and 183 GeV, taking into account common systematic errors in the 
theoretical prediction and in the event selection. The theoretical prediction of o"ww is 
computed with GENTLE ||18|| , taking into account its residual sensitivity to a change in 



the total width. The possibility that the triple gauge boson couplings could be anomalous 
is allowed, by fitting simultaneously Ak 7 , A 7 and Agf within the constraints allowed by 



the analysis of angular distributions in the semileptonic channel ||22|| . The W mass is 
allowed to vary in the fit within the constraint imposed by the current world average 
value (not including the W mass from the WW cross section measurements at LEP), 
m w = 80.43 ± 0.08 GeV/c 2 J[9|. 



In addition to the systematic errors on the cross section measurements, three more 
correlated systematic errors are considered: i) a theoretical uncertainty of ±2% in the 
prediction of the total WW cross section, ii) a theoretical uncertainty of ±0.5% in the 
Standard Model W width and iii) the beam energy uncertainty AE^ = ±25 MeV, where 
the last two sources of systematic uncertainty have been found to have a negligible effect 
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(below 1 MeV). The different contributions to the total systematic error are summarized 
in Table [3[ 

The fit to the cross section data jlj gives 

I™ = 30 til (stat.) ± 33 (syst.) MeV (x 2 /dof = 3.5/2), 

T^ vis < 139 MeV at95%C.L., 
B invis < Q 5% at g5% c L 

The corresponding total W width is 

T w = 2.126 IK (stat.) ± 0.035 (syst.) GeV. 



Table 3: Summary of systematic uncertainties on the T^ V1S determination from the cross 
section measurement. 



Source 


AI™ [MeV] 


Cross section measurement 
Theor. uncertainty on cr\vw 


±24 
±23 


Total uncertainty 


±33 



7 Conclusions 

Using an integrated luminosity of 57 pb" 1 the W pair production cross section at ^fs 
= 182.7 GeV has been measured in all decay channels. The determination of individual 
branching ratios has been performed. The total cross section is 15.57 ± 0.62(stat.) ± 
0.29(syst.) pb. After adding the data taken at 161 and 172 GeV cm. energy, the 
hadronic decay branching ratio is 68.93 ± 1.21(stat.) ± 0.51(syst.)% which is used to 
determine the CKM matrix element \V CS \ = 1.043 ± 0.058(stat.) ± 0.026(syst.). These 
results are consistent with similar measurements performed at LEP at the same centre- 



of-mass energy p3| . 

The agreement of the cross section with the Standard Model prediction allows a limit 
to be set on the W partial width into invisible decays: T^ V1S = 30^g (stat.) ±33 (syst.) MeV, 
r™ vis < 139 MeV at 95% C.L. 
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Appendix A. Hadronic neural network input variables 

The neural network hadronic event selection uses 19 variables. These are based on global 
event properties, heavy quark flavour tagging, jet properties and WW kinematics and 
are listed below together with their discriminating power. The four jets are numbered in 
order of decreasing energy. 
Global event properties 

• Fox Wolfram moment HO (4.4%) 

• Fox Wolfram moment H2 (4.7%) 

• Fox Wolfram moment H4 (10.1%) 

• Sphericity (3.9%) 

• Missing Energy (4.4%) 

• Sum of the pf of all charged tracks to the beam axis (4.1%) 
Heavy flavour tagging 

• Sum of the b-tag probabilities for the four jets (5.7%) 

The following jet related variables are determined from kinematically fitted jet 
momenta. 

Jet properties 

• Number of good charged tracks in jetl (6.2%) 

• Maximum energy carried by one energy flow particle in Jetl (3.8%) 

• Maximum energy carried by one energy flow particle in Jet2 (4.6%) 

• Maximum energy carried by one energy flow particle in Jet3 (4.7%) 

• Sum of the angles between the leading and remaining charged tracks in jetl (5.5%) 

• Sum of the angles between the leading and remaining charged tracks in jet2 (3.6%) 
WW kinematics 

• Sum of the cosines of the 6 angles between jets (8.7%) 

• Cosine of the angle between jet2 and jet3 (4.6%) 

• Energy of jetl (8.1%) 

• Energy of jet2 (3.8%) 

• Momentum of jet4 (5.5%) 

• Asymmetry in momentum between jet2 and jet3 (3.8%) 
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